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ABSTRACT 


The results of the measurements of carbon dioxide in 
air at Barrow, Alaska, are presented, Reference gas ccom- 
parison data are tabuleted, and methods of calculations 
are discussed. The average daily concentrations of 
atmospheric carbon dioxide are tabulated for the period 
10 July 1961 to 15 September 1967. The diurnal variations 


of carbon dioxide during this period are also presented. 
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PREFACE 


This ‘eport presents a summary of measurements of the concentration 
of carbon dioxide near Barrow, Alaska, and at the North Meadow Lake field 
station of the Arctic Research Laboratory, Barrow, Alaska. 

This work was supported under a contract from the Office of Naval 


Research (ONR 477(24)) with the Department of Atmospheric Sciences, 


University of Washington. The analytical program was conducted cooperatively 


with Dr. C. D. Keeling's (The Scripps Institution of Oceanography, La 
Jolla, California) atmospheric carben dicxiic program at Mauna Loa, 
Hawaii, and the South Pele Station, Antavetiva. 

The carton dioxide trogrmi at Bacrow win Snitiated and maintained 
from July 1961 to fugust 2967 by J. J. Keliny, Jr. Operation of the carbon 
Gioxide ausiyzer wan soetineed by dU. Unger fron September 1962 to March 
W9, ang from April Ges ty Gctener 1960 by o. Stroschein. 

Tee carbon dioxide progras at tha korth Headow Lake field station 
Yas oof atcsinec hv f, weaver frost Januvary L965 to September 1965, by 
% Sates fey September 1G66 ts Septembes 2G66, and by W. fowald from 
O.toter }#bb to the conclusion of the crograsy iv September 1967, 

The Director af tae Awlic fmaearch Lekoratory, Dr. M. C. Brewer 


gad his staff, have prowidec invaluaile asucstance to the project. 
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I. INTRODUCTION 


This report presents the final results of a program to measure the con- 
centration of carbon dioxide in the atmosphere at Barrow, Alaska from 
July 1961 to October 1963, and at the North Meadow Lake field station 
of te Arctic Research Laboratory, Barrow, Alaska from January 1965 
to September 1967. The location of each site is shown in Figure l. 

Dailv average values of the concentration of CO in the surface air 
are computed from da.a derived from original strip chart records of a 
continucus recording infrared gas analyzer installed at the station. 
All relevant data and computations for the vears 19641 through 1967 
are contained in this report. The experimental method is described 
in a Technical Report by Kelley (1964). Interpretation of the data 
will be published in a scientific journal. 

All] data in this revort are final values. The procedure for 
computation follows that used at the Scripps Institution of Oceanography, 


La Jolla, California as described in Research Renorts I throuph VITI*. 


If. REFERENCE CAS COMPARISONS - TABLE J, la 


Yhis table lists, in chronological order, beth the observed scale 
differences which were used to calibrate the eas analvzer, and their 
conversion into index valucs proportional to fA, concentration. The 


kConies of these reports may be obtained from Pr. Charles PD. Keelinr, 


§.I1.0., La Jolla, California. 
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calibrations of the analyzer consisted of repeated comparisons of pairs 
of apecially prepared gas mixtures of CO, in nitrogen obtained from the 
Scripps Institution of Oceanography. These reference gae mixtures were 
stored in stainless steel cylinders called "canka." 

Under normal operating conditions ten comparisons were obtained by 
alternately passing one gas of the pair, and then the other, through 
the infrared analyzer for five minutes at the same flow rate employed 
in the air measurements (normally 0.5 liters per minute). As soon as 
one series of tcn comparisons was run, one or both tanks were replaced 
and another pair of tanks compared. This process was repeated, as a 
rule, three times during the life of a “working reference" tank: at 
the beginning of use, when the gas pressure was half depleted, and at 
400 p.s.i., before return of the cylinder to Scrippe for final calibration. 
The scale difference between two successive traces, in recorder chart 
ordinates, was read with a straightedge scale (30 divisions to the inch) 
by drawing parallel straight lines through each of the traces belonging 
to the reference gases. The successior of individual scale differences 
for each tank pair was entered on data sheets. The calibrations described 
served to estabiish the recorder sensitivity of the infrared gas analyzer. 
They were also used to determine the index values of reference gases 


known as “working references," used in connection with air mesecurementa. 
The general method used by the Scripps Institution for tank standardization 
(Figure 2) was modified somewhat ai Barrow during the period covered 

by the report. Initially, two reference gases were compared with air 


every 30 minutes. Frou 1962 to the end of the program in 1967, three 


reference gases were employed. Two of the reference gases were used 
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as comparison standards. The third reference gas was ranked as a 
“working reference" and compared with air every 30 minutes. The "working 
reference" was compared directly to the comparison standards, primary 

and span (high or low), in the absence of a suitable secondary reference 
at Barrow. 

In Table 1, tanks related to each calibration run are identified by 
tank numbers. They are listed in columns 1 and 2. Column 1 lists the 
number of the standard tank; column 2 lists the number of the tank with 
which it is compared. The average observed scale difference for each 
tank pair is listed in column 3. A positive number indicates that the 
compared tank, for that particular comparison, has a higher scale reading, 
and consequently a higher CO, concentration than the standard; a negative 


number indicates the reverse. The number of comparisons (not always 


10) which entered into each average scale difference is listed in column 4. 


The reccrder sensitivity, determined by comparisons of the primary 
and secondary standards, is expressed by a "recorder scale factor" listed 
in column 5. The computed index differences, each with the same sipn 
as the corresponding scale difference, and the index value of the compared 
tank, in terms of a prescribed, or "assigned" value of the standard tank, 
are listed in columns 6 and 7. 

Weighted average recorder scale factors, column 5, are copied from 
colum 7 of table 2. Index differences, column 6, were computed by * +e 


formula given in Section III-D. 


So Aeaesietatinneath 


The computed index values, colum 7, are the algebraic sums of the 
index values of the standard tanka listed in colum 1 and the index 
differences of column 6, The index values of the standard tanks are those 
listed in column 7 of Table 6. In the case of the primary tanks designated 
as I and Ipg or Ipc in Figure 2, they depend solely on the measurements 
made at Scripps. 

The data listed in Table 1 are for the period 2? January 1965 to 
18 September 1967 at the North Meadow Lake field station. The data in 
Table la are for the period ] April 196% to 14 June 1967 at the University 
of Washington. The data for reference gas comparisons at the Arctic 
Qcean beach site from 6 January 1962 to 21 August 1963 are given in 
Report 2 (Kelley, 1966). 

The scale differences entered in column 3 of Table 1 and la are 


copied from original entries on the reference gas data sheets. 
III.RECORDER SCALE FACTORS - TABLE 2, 2a 


A. Definition 

The Recorder Scale Factor, RSF, is defined as the index difference 
between two reference gases divided by the number of scale divisions 
(30 div. to the inch) between the recorder chart traces for each reference 


gas recorded in units. Ten divisions on the scale were equal to one unit. 


The index of a reference gas tank is defined as a provisional C0, 


concentration, or mixing ratio, in ppm based on the initial analysis of the 
tank. 
B, Standard Computation for Three Mutually Compared Tanks 

Columns 1 through 4, except values with asterisks explained below, and 


the last entry in column 4 for each calibration day, list selected data 


2 thie 


copied directly from the corresponding colwms of Table 1. These data are 
employed, as shown below, to obtain daily RSF's. 

The following format has been adapted from the Scripps (SIO Reports 
1-8) scheme for tank standardizations, where the standard tank index 
values have been assigned the symbcls A and H, and the compared tank, 


without such assignment, the symbol X: 


Standard Compared Observed Scale No. of 
Tank No. Tank No. Difference Comparisons 
A X (X] - [A] a 
B x (x) - (8) b 
A B (y)}* (a or b)* 


[Y} = C UX} - (A}) - C(X] - (8) = C (By - (A) 


The asterisk on (YJindicates the calculated value, and the brackets 
indicate index values. The number of comparisons assigned to [Y) is "a" 
or "pb", whichever is smaller. The observed scale differences were taken 
from the averages of the individual chart scale differences (Table 1). 

The observed comparisons and calculated values of any group of three 


tanks, A, B, and X are set off in the table by boxes. 


C. Determination of Index Differences 
Index differences shown in Tables 2 and 2a, column 5, were obtained 


from index values of the standard tanks as follows: 


Tank Numbers Index Values Difference 
11589 vs 18206 314.59 - 310.70 3.89 
11589 vs 18207 314.59 - 312.74 1.81 
18206 vs 18208 314.59 - 312.40 2.19 
11£89 vs 10071 -314.59 + 339.00 24.41 


a Bt 


The index values of separate tanks are obtained from Table 4. 
D, Weighted Average Recorder Scale Factors 

Column 6 presents RSF's computed according to the definition given 
in section III-A. Column ? presenta the weighted average values of the 
recorder scale factors for each calibration day. ‘The index differences 


in column 6 of Table 1 are computed by the formula: 


Computed Index = (Observed Scale Difference) (Weighted Average RSF) 
Difference 


IV, SUMMARY OF RECORDER SCALE FACTORS - TABLES 3, 3a, 3b, 3c 


en eR re mn 


Values of RSF's are assembled in chronological order in column 
5 of Table 3 for Periods 1 through 11, 10 July 1961 to 16 December 1961. 

A graph for the data from Period 1 through Period 11 is shown in Figure 3. 
During these periods, breaks in continuity were extreme due to the 
development of a method for determining a standard operating method for 
infrared analysis at the Arctic Ocean beach site. Each break is charac- 
terized by a change in the span control on the infrared analyzer. There- 
fore, discrete periods of relatively constant recorder scale factors are 
created between span changes. 

During periods 12 to 15, from 6 January 1942 to 21 August 1563, 
(Table 3a) a constant recorder span setting was used and a mutual comparis«r 
method of tank standardization was set up. These data are snown in 
Figures 4 and 5. 

Table 3b presents a summary of recorder scale factors for the period 


2 January 1965 to 18 September 1967 at the North Meadow Lake field station. 
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The data are derived by means of the mutual comparison method. Discon- 
tinuities in the recorder scale factors oceurred at various times during the 
period of analysis at North Meadow Lake. Table 3c provides daily 
recorder scale factors for the periods where discontinuities occur. 
The data in Table 3c were produced by drawing a line between the end 
point of the average RSF curve for one period to the starting point of 
the average RSF curve of the next period. Daily RSF's were recorded and 
defined as eliding recorder scale factors. The recorder scale factors 
adjusted for pressure barometric for the period 2 January 1965 to 18 
September 1967 are shown in Figures 6, 7, and 8. 

All recorder scale factors were adjusted to a standard barometric 
pressure of 30.00 inches of mercury (1015.9 mb). The following formula 


was used: 


Adjusted Recorder, (Weighted Average RSF) (Observed Barometric Pressure) 


Scale Factor 30.00 inches of mercury 
V. INDEX VALUES OF REFERENCE GASES - TABLE 4 


This table presents the index values of all reference gases used 


emma ne i tt AA I AeA AEE Etat 


in the Arctic field program and at the University of Washington. All 


of the analyses were performed at the Scripps Institution of Oceanography. 


VI. COMBINED SCRIPP! AND BARROW INDEX VALUES 
OF WORKING REFERUNCE GASES - TABLE 5 
This table sunmarizes the index values of all working reference 
gases used at the North Meadow Lake field station from 1965 to 1967. 


The data for working reference gases for 1961 to 1963 are given in Report 


2, Table 7 (Kelley, 1966). 
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Entries in columns 2 and 3 are the numbers of comparisons and weirhted 
averepe index values based on measurements at Scripps prior to and after 
use. These data are taken from Table 4. 

Trertes in columns 4 and 5 are the numbers of comparisons and weirhted 
averare index values hased on measurements at the North ‘eadow Lake 
field st.tion near Barrow, Alaska. 

Fntries in columns 6 and 7 list the total number of connarisons and 
veiphted averape index values hased on all measurerents at Scripps and 
at TRarrow, 


VII. COMPARTSON OF SCRIPPS AND RARPOW INDFS 
VALUES OF WORKING RFFIPINCE CASTS = TALI &. 


This tahle presents the results of all analvsis of Parrow reference 
tanks used in the eaily comnartson with afr during 195 ta 2947, This 
comnarison for tne vears 194} to 1°43 is riven in “Nenort 2, Tahic & 
(Yellev, 1964), The table comnares the results of the analyses made 
at the Scrinns Institution of Oceanopranhy with those rade at Marrow, 

The total number cf indices and comnarisens shown {n columns 2? and 3 are 
tai:en from the columns headed cumulative runs and cumulative index of 
Table 4, The numher of cornarisons and reference tank indices in coJumns 
4 ane 5 are conied from columns 4 and 5 of Table §&. Intries in column 

6 are the differences hetween index values in narts ner mi}lion he volume 
CN, ohtainec from measurements at PRarrov and Scrinns, and represent 


the {tndey denarture from Scripps values. These cata are nletted in Fieure 9 


for tire neriod 1°6)-1963, and Firure 10 for the nertod J°45-1967, 
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VIII. INDEX VALUES OF AIR - TABLE 7, 7a 


A. Computation of Index Values 

These tables contain all daily co, values at the Arctic Ocean 
beach site and the North Meadow Lake field station for the period 
1961-1967. Colwnn 1 lists the calendar date of the analysis. Column 2 
shows the observed daily average scale difference between the air trace 
and the working reference gas trace. Colum 3 lists the total number 
of comparisons for the day. A full day consists of forty-eight half- 
hour comparisons. Column 4 records the average barometric pressure for 
the day. These values are obtained from the original carbon dioxide 
daily data sheets. All of the data for the years 1961 to 1963, Table 7, 
were processed manually. All of the data from 1965 to the termination 
of the program in 1967 were precessed with the aid of an IBM 7094 
computer. The Fortran 4 format is given in Appendix 2, 

The observed scale difference in column 2 was adjusted to a standard 
barometric pressure of 30.00 inches of mercury (1015.9 mb) by the 


formula: 


Adjusted Scale . Observed Scale , 30.00 inches of Hg 
Difference Difference Average pressure for the day 
and this adjusted scale difference is recorded in colum 5. 
Column 6 lists the RSF's to be used in converting scale uifferences 
to index units. These are based on values listed in Tables 3, 3a, 3b, 
and 3c which have been adjusted to standard atmospheric pressure. 


The computed index difference in colurmm 7 is obtained by: 


Computed Index Difference = RSF X Adjusted Scale Difference 


- 10- 


The working reference tank used in the comparieon with air is listed 
in colum 8, its index in colum 9. These data are taken from Table 5. 
The daily average index values in column 10 are the algebraic sume of the 
entries in columns 7 and 9, 

In several instances reference tank changes were made during the 


day. In these cases the daily index was computed by first computing 


an air index for all of the observations with one tank in use, and doing 
the same for the data taken with the replacement tank. These data are 
reported separately in column 10. The daily index is taken as the average 


of the two separate index values. 


B. Manometric Concentration Scale 

Throughout this report, the co, data have been reported in terms of 
an index scale. It was established provisionally at Scripps in 1959 
that the true concentration in parte per million by volume is related to 


the isdex scale by: 


Manonetric Concentration = (C - 311.51) 1.2186 + 311.51 


This equation is based on the absolute calibration of primary and 
span reference gases by means of a mercury manometer. It is possible 
that future manometric calibrations may result in further adjustment of 


the index scale. The intercept value (311.51) is believed to be correct 


within 1 ppm; the slope value (1.2186) within 0.03. 


This equation has been used to report atmospheric co, concentrations 
in the reports and journal contributions listed in Appendix 1. 

Column 11 lists the manometric concentration of the index values 
reported in column 10 in ppm co, by volume dry air. The data are shown 
graphically in Figure 11 for the period 1961-1963, and Figure 12 for 


the period 1965-1967. 


1X. MONTHLY AVERAGE INDEX OF 
CARBON DIOXIDE - TABLE 8, 8a __ 


This table presents the monthly average index of atmospheric CO, 
at Barrow, Alaska for the years 1961-1963 (Table 8), and 1965-1967 
(Table 8a). Columns 2, 5, and 8 list the number of days for which air 
index values are quoted in Tables 7 and 7a. Columns 3, 6, and 9 list 
the monthly carbon dioxide index values. 
X. MONTHLY AVERAGE INDEX OF CARBON DIOXIDE, 
MANOMETRIC CONCENTRATION SCALE - TABLE 9, 9a 
The data in Tables 9 and 9a are the manometric indices in ppm 
co, for index values quoted in Tables 8 and 8a. 
The averages of the monthly values for each year are given in 
columns 3, 6, and 9, 
XI. TWELVE MONTH RUNNING ME/.N CONCENTRATION 
OF CARBON DIOXIDE - TABLE 10 
A twelve-month running mean concentration of CO was calculated to 
smooth out the seasonal variation in CC). Concentrations for co. 
are presented for the period 1962 to 1967. Means are plotted in Figure 
13 versus the sixth month of the appropriate 12-month interval. Least 


squares regression lines were fitted te the data for the periods 1962- 


1963 end 1965-1967. Equations for each regression line were determined as 


-12- 


cn Pn nn ret TO I ET EEE OE 


2 - 


o.oo ee a ee eget 
eos. ve ca é 1962 
ae ve . : 2 SANS 4 
oe a OS x sel 
2 | 7 “- “A 
je) i. $ 
pa Ie \y ‘ i e 
a : : 
Baa ae 
= 310 iz in i o . 4 
uJ : acme a | 
= oe : 
oO 1 
Cp Oe 1 
w g ; 
Ce oes ee aie hs, Rd ls ed 
JAN FEB MAR APR MAY JUN JUL AUG SEP ocT NOV DEC 
375 _— i 7 | 7 
1 1 T ! { i | 1 ‘ 
: +, a 
ook, jot 968 | 
0K. ty . 
Se ee My: i 
| . ey i 
\ Be : 
315 7. ; 7 — 
a . 
310 L nie © ae 7 
SS 
305 
| 
cred eens Me! ieee SE GO: SE ae tee 
JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC 
Figure 11 Daily Average Concentration of co, 1961-1963. 


325 eas Sa ae es 


i961 
320 }- : 

; Lig 
eI ; ee oes 

r e . : 
310 LL me ‘ ve 4 

Slee | 

] . aN a | 
305 +- 1 

| 


Fr) Lies Caney Coe LOR SSE ire? ee Boe 
JAN FEB MAR APR MAY JUN JUL AUG SEP OCT ie OEC 


oe ee ee CRs a a 

. Re os < 1965 3 

lind yay %: l 

20 EE : ~j 

o F < of REO 

Se cae | 

a1 a 5 
i ene 

310 t- Se 7 

| 

1 

305 E a 

| | 

H | 

fe Sa Gop ths Spee fhe ee he ly 


JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC 


‘ee Pert Tied oe ge 
.. sb meee PY 
es ~~ e a ah oy 

~~ 320 as * P> NIE 

! 8 ‘ eee ee | 

z | * more | 

O ws. 

ae yN5 — op me 7 

a i net : a: 

a | *- = . ‘ 

Fe Cas : 

= sOr- one a 

uJ 

oO | | 

z 

S ws '- a 

7 | 
iA 

yp a Th 

JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC 
eC ia Le ea oD ds eae ed aod 
a ae eS 


300 st tp eg Ef ae + ee | J 
JAN FEB MAR 4PR MAY JUN JUL AUG SEP OCT NOV CEC 


Figure 12 Da. ly Average Concentration of CO 


2 1965-1967. 
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Comusite averages appear in column 5. Column 6 lists the departure 
of these sverages from the annual mea value of 314.69. The monthly 


average concentration of C0) near Barrow, Alaska referred to a constant 


datum (1960) is shown in Figure 14. 
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XUIT. THE DIURNAL VARIATION OF CO - TABLE 12, l2a, 13¢ 


The diurnal variation of co, was calculated from the half~hourlv 
observations recorded on the oririnal carbon dioxide data sheets. 
The dturnal variations are based on hourly air indices by finding the 
averape index for each hour of the dav for each month. 


Firure 15 and Table 12¢ show the averape diurnal variation of 


CQ for each month hased on the data for CO) from 1961 to 1967. 
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Figure 15 Average Diurnal Variation of co. Near Barrow, Alaska, 
1961-1967. 
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APPENDIX 2 


FORTRAN 4 COMPUTER PROGRAM 


This program is designed for use cn the IBM 7094 computer. Data 


for ca, from five sources of air are entered into the program for each 


index for each half hour interval is 


is given for each day as follows: 


Number of haif hourly observaticns 


Recorder Scale Factor 


Mean Index 
Mean Index 
Mean Index 
Mean Index 
Mean Index 


Mean Index 


1 


2 


3 


4 


5 


all levels 


Mean barometric pressure 


> 


Mean temperature, ~C 


Location of leveis 1 through 5 


INSTRUCTIONS FOR KEYPUNCHING DATA 


1. The first data card 


the data of the day's run, e. 


has the form: three blank spaces, then 
ge: 


17 _ March _ 1965 


Otner data may also be placed on this card and will be reproduced in the 


i 

t 
half hour interval. The co, 
computed. Wind direction, wind speed, barometric pressure, and ambient 
air tempe‘ature are logged in the data output. A summary of the data 
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Daseably 


heading of the output just as printed on the car. 

2. The second card contains the chart number (3 spaces), reference 
tank number (6 spaces), and the index (8 spaces) 

For example, if chart number 12, reference tank, 11633, and the 


tank i.udex, 310.68, are the data, then the card reads as follows: 
0120116330310.680 


3. Data Cards: There must be a total of forty-eight cards; 
i.e., one for each half hour of data even if data are missing for part 
of the time. The pressure, temperature, wind direction, and wind 
speed data must start on the first card. All data must start on the half- 


hour and alternate from there. For example: 


Wind Wind 
Time Ll L2 L3 Lu LS Press. Temp. Dir. Speed 
0030 2.93 7.82 10.63 15.72 2.93 130 13 1,029 0 


0100 2.83 7.89 11.02 15.47 2.83 - - - - 
Cards read: 


00?.83907.89011.01015.47002.8399.39 


4, Missing Data: 

a) If more than one cf the index values are missing, punch 99.99 
in the first five columns. Then skip to column 31 and punch in the wind, 
temperature and pressure data on the half hour. 


b) If only one of the index values is missing, substitute an average 


value for the missing index value. 


c) If the pressure, temperature, and wind observations are missing 
for a particular half hour, punch 99.99 in columns 31-35. 
S. Program Format: 


Tne Fortran-4 program and a specimen of the arrangement of the 


data for one day are given. 


SJOB 3,5,6000 2002650 
SFORMS OUL PROGRAM , 141326 
SEXECUTE I3JOB 
SIBJOB GO ,MAP 
SIBFTC COTW1 
c CARBON DIOXIDE INDEX PROGRAM-SUMMER 5 LEVEL VERSION(REVISED) 
DIMENSION Si(S,48),AI( 5,48) ,KT(48) ,SDAC 48) ,T624),P( 24) IWD( 24), IWS 
1(24) 
l FORMAT(13,16,F8.3) 
2 FORMAT(SF6 .2,F5.2,F4.0,213) 
10 = FORMAT(1H1,40X,38H DAILY CARBON DIOXIDE AIR INDEX VALUES/53X,15H B 
LARROW, ALASKA/&@¥,26H NORTH MEADOW LAKE STATION//) 
11  FORMAT(89H 
1 ) 
12. FORMAT(19X,15H CHART NUMBER 213,15%,18H REFERENCE GAS W °16,8X,8H 
INDEX 2F9.3/) 
13. FORMAT(6H AST,10X,245H HALF HOURLY SCALE DIFF.,19X,23H HALF HOUR 
1Y INDEX, AIR,13X,35H WIND OIR SPEED PRESSURE TEMP/10X,3H 1., 
24X,3H 2.,4X,3H 3.,4X,3H 4.,4X,3H £.,6X,39H 1.,7X,3H 2.,7X,3H 3.,7%, 
33H 4.,7X,3H 5.,4X,33H DEGREES KN. BARS c.) 
14 FORMAT (2X ,24%,1X,F7.2,F7.2,F7.2,F7.2 577+ 2,4X F8.39, 2% FB .3,2% F869 52 
1X ,F8.3,2X,F8.3,5X,13,6X,19XUX ,F7.4,3X,F5.0) 
1S FORMAT(2H I4&,16H DATA MISSING) 
16 FORMAT(//37HO NUMBER OF HALF HOURLY OBSERVATIONS=F4.0/25HO RECORDE 
1R SCALE FACTOR 2F7.3/16HO MEAN INDEX 1 = F8.3/16HO MEAN INDEX 2 =F8 
2.3/16HO MEAN INDEX 3 = F8.3/16HO MEAN INDEX 4& = F8.3/16HO MEAN INDEX 
3.5 = F8.3/28HO MEAN INDEX AT ALL LEVELS = ¥8.3/17HO MEAN PRESSURE 27 
47.1,10H MILLIBARS/20HO MEAN TEMPERATURE =F5.0,19H DEGREES CENTIGRA 
5DE/16HO L1=16.0 METERS/15HO L2#1.0 METERS/15KO L3=0.5 METERS/16HO 
6L420.2S METERS/17HO L5=0.125 METERS) 
17 FORMAT(2H = =14,16H DATA MISSING, 77X,13,6X,13,4X,F7.4,3X,F5.0) 
18 FORMAT(1H1,114X,7H PAGE 2) 
RSF=1.720 
Kl=¢ 
DO 50 J=1,48 
MI=J 
IF(MOD(MJ ,2))51,52,51 
51  K1=K1+30 
KT(J )#Ki 
K1=K1+70 
GO TO 50 
52. KT(J)=K1 
50 CONTINUE 
200 ITR=0 
K1=0 
K2=0 
READ(S,11) 
READ(S5,1) ICN,IRT,WRG 
DO 75 J=1,48 
Kizv+l1 
J1lsJ 
IF(MOD(J1,2))76,77,76 
76 K2=K1-ITR 
READ(S,2) (SD(I,J),I=1,5),P(K2),T(K2),IWD(K2) ,IWS(K2) 
ITREITR+1 


GO TO 75 


AV320.0 
AV4&=0.0 
AVS=0;:0 
; DIN20.0 
| DO 100 J=),48 
A=SD(1,J) 
IF(A-99.99)i01,100,101 
103 DIN#DIN+1.0 
AV1Z£AV14SD(1,9 ) 
AV 22AV24+SD(2,0 ) 
AV3=AV3+SD(3,J) 
AVUZAVU+SD(4 0) 
AVS=AVS4SD(5,J ) 
| 100 CONTINUE 
AV1=AV1/DIN 
AV2=AV2/DIN 
AV3=AV3/DIN 
AVUSAV4U/DIN 
AVSsAVS/DIN 
AV=0.0 
DO 102 J=1,48 
A=SD(1,J ) 
IF(A-99.99)103 ,102,103 
103 S$=0.0 
pO 104 I=1,5 
AI(I,J )=WRG+(SD(1,J )*RSF ) 
104 S=S4SD(I ,J) 
SDA(J)=S/5.0 
AV=AV+SDA(J) 


j 

77 READ(S,2) (SD(I,J),I=1,5) 
75 CONTINUE 

DO 99 Jul,48 

SDA(J)#0.0 

pO 99 I=1,5 
99 AI(I,J)20.0 

AV1=0.0 

AV2=0.0 


102 CONTINUE 
AV=AV/DIN 

| AV1=9RG+ (AVIFRSF ) 
AV 2=WRG+ (AV24RSF ) 
AV3=¥RG+(AV3"RSF ) 


| AV&=WRG+ (AVUSRSF ) 
AVS=WRG+(AV5®RSF ) 

AV=WRG+ (AVARSF ) 

DO 105 Y=1,24 
| T(1)=0.5568(T(1I)-32.9) 

IF(P(I)-15.00)316.316 ,317 
316 P(1)=P(1)/10.) | 
GO TO 105 
317s FC I) =(33.864*P(I))/1000.0 
105 CONTINUE 

PA=0.0 

TA=0.0 

DIT=0.0 

DO 125 I=1,24 


Te eae = a 


A=P(1) 
IF(A-99.99)226,125,126 

126 DIT2DIT+1.0 
PAsPAtP(1) 

TA®TA4T (I) 

125 CONTINUE 
PAx(PA*®1000.)/DIT 
TASTA/DIT 
WRITE(6 ,10) 

WRITE(6 ,11) 

WRITE(6,12) ICN,IRT,WRG 
WRITE(6 ,13) 

Kiz0 

DO 150 J21,48 

Ax8D(1,J) 
TF(A-99.99)151,152,151 

152 JlsJ 
IF(MOD(J2,2))155,156,155 

155 K12K1¢1 
B=P(K1) 
IF(B-99.99)157,156,157 


157 WRITE(6,17) KT(J),IWD(K1),IWS(KL) ,P(K1) ,T(K1) 


GO TO 1580 

155 WRITE(6,15) KT(J) 
GO TO 150 

151 J1eJ 


IF(MOD(J1,2))153,154,153 
153 Kl=Kit+l 


WRITE(6,14) KTiJ),SD(1,J),SD(2,J),SD(3,7),SD(4,7),SD(5,J),AT(1,J), 
1AI(2,0) ,AI(3 50) ,ATC4,5) AICS 0), IWDCK1) , TWS(KL) ,P(K2) ,T(KL) 


GO TO 150 


154 WRITE(6,14) KT(J),SD(1,7),SD(2,0),SD(3,7),S0(4,5),SD(5,J) AI(1 J), 


1A1£2,5),AT(3,5),AICY J) ,AT(5 J) 
150 CONTINUE 

WRITE(6,18} , 

WRITE(E 22 


WRITE(6,16) DIN,RSF,AV1,AV2,AV3,AV4,AVS ,AV,PA,TA 


GO TO 200 
END 
SENTRY COTW 
SDATA 
1 MAY 1967 
1.280100670311.71C¢ 
002 .93007 .82019.63015 .72002.9310.29032.130013 
0 2.83007 .89011.01015.47002.8399.99 
002 .96007 .75011. 26015. 30002.9610.29032.140011 
002. 93007 .50011.05015.18092.9399.99 
003 .CB007 .61009 .02015.72003.0810.29037.140011 
002. 90007 .75010. 96016 .17002.9099.99 
002.69006.09010.75016.76002. 8910. 29031].156010 
0 2.89098. 20010 .87017.39002.8999.99 
002 .83008.57008.57017.88002.8310.28031.14%0011 
002.93008 . 87012 . 23017 .93002.9399.99 
002.94908.91011. 35018. 29002.9410.26031.14U009 
002 .88009.21021.13018.54002. 8899.99 
002.90008.66011. 21019.23002.9010.28032.2140008 
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